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Fluctuation properties of the Langevin equation including a multiplicative, power-law noise and a 
quadratic potential are discussed. The noise has the Levy stable distribution. If this distribution is 
truncated, the covariance can be derived in the limit of large time; it falls exponentially. Covariance 
, in the stable case, studied for the Cauchy distribution, exhibits a weakly stretched exponential shape 

^— ( ' and can be approximated by the simple exponential. The dependence of that function on system 

' parameters is determined. Then we consider a dynamics which involves the above process and 

, obey the generalised Langevin equation, the same as for Gaussian case. The resulting distributions 

» I ■ possess power-law tails - which fall similarly to those for the driving noise - whereas central parts 

' can assume the Gaussian shape. Moreover, a process with the covariance 1/t at large time is 

, constructed and the corresponding dynamical equation is solved. Diffusion properties of systems for 

. both covariances are discussed. 

in ; 

I \'. I. INTRODUCTION 

^ . 

D I Trajectories encountered in complex systems often reveal discontinuities and the probability distributions are not 
' governed by equations with local operators. Those distributions differ from the Gaussian and contain slowly falling 
I , tails. Levy stable distributions with tails ~ |x|~^~", where a € (0, 2), are distinguished due to the generalised central 
"j^ ' limit theorem. However, the divergent variance makes the Levy stable distributions problematic in some physical 
applications; it may imply, for example, infinite kinetic energy. Similarly, covariance functions for the Levy stable 
processes with a < 2 do not exist. The above difficulties do not emerge if the tails, being still of the power form, fall 
I faster than for the Levy stable distributions. In fact, tails of the form |a;|~'^, where /3 > 3, are frequently observed. 
' This is the case for the financial market that possesses typical characteristics of the complex system and then some 
I of its properties may be universal. Analyses of returns of stock indices show that a cumulative distribution of returns 
' O . is power-law with 3 < /3 < 4 [l|, 01 and such values of /3 are required by the optimal market strategy Q . The minority 
^ ■ game implies a similar value, /3 = 3.9 Since large jumps represent extreme events, one can expect that first 
passage time probability should obey the WeibuU distribution. However, a phenomenological analysis of the empirical 
data demonstrate that it is the case only for small returns, the large ones are of the power form with (3 = 3.32 [|. 
' Moreover, fast falling power-law tails result from a multifractal analysis of the extreme events JiJ, characterize the 
I hydraulic conductivity in the porous media and the atmospheric turbulence Q- 

' The variance becomes finite when we modify the asymptotics of the Levy stable distribution by introducing either 
. a simple cut-off or some fast-falling tail. Such truncated distributions very slowly iS] converge with time to the 
' normal distribution. Moreover, dynamical systems stimulated by the Levy stable noise possess finite moments of the 
fS| , stationary distribution if the particle is trapped inside a potential well with a sufficiently large slope Q . The Langevin 
■ equation with a multiplicative Levy noise t] and a linear deterministic force also predicts finite moments, if interpreted 
I in the Stratonovich sense; this property was demonstrated by numerical simulations [l^ |Tl[ . On the other hand, we 
I can define the coloured noise r]c{t) by the generalised Ornstein-Uhlenbeck process, 

T1 ■ dr]c{t) ^ -jnTlc{t)dt + jnd.L{t), (1) 

where the increments of L{t) have the stable Levy distribution, and take the white-noise limit, ri(t) = lim^^_>oo Vcit)- 
Stochastic equation with the white noise rjjt), given by the above expression, allows us to change the variable in the 
usual way and obtain the Stratonovich result I'M ■ Tfie generalised Wiener process with that noise as a driving force is 
characterised by the subdiffusive motion. Taking into account that the mass in the Langevin equation is finite modifies 
the slope of the tails: it diminishes with the inertia and finally converges to the result of the Ito interpretation for 
the infinite mass. Multiplicative non-Gaussian white noises serve to describe population and ecological problems: a 
dynamics of two competing species [ij, and a population density in terms of the Verhulst model Jil5 Jl6|. 
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Long tails of the distributions in complex systems are often accompanied by a long memory [J, [17[ and then 
time-dependence of the fluctuations becomes non-trivial. The well-known multifractal structure of the financial time 
series J^] is attributed to both fat tails, that fall as a power-law but faster than for the stable Levy, and power- 
law correlations [l^. A slow decay of the correlations, even slower than a power- law, is regarded as a necessary 
condition of the multifractality in any complex system [20j. In the present paper, we demonstrate that the generalised 
Ornstein-Uhlenbeck process £,{t), driven by r]{t), possesses a well-determined covariance and fat tails of the distribution. 
Moreover, trajectories ^(i) preserve a typical feature of the Levy flights: smooth segments interrupted by large, rare 
jumps. 
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In the presence of the memory effects, a description in terms of the standard Langevin equation with a coloured noise 
is problematic since a response of the system to the stochastic stimulation is not instantaneous and, as a consequence, 
a retarded friction must be introduced. The fluctuation-dissipation relation requires that the equipartition energy 
rule is satisfied, i.e. the temperature is well defined, and a retarded friction kernel is uniquely determined by the 
noise covariance; then £^{t) is called the internal noise. Obviously, that relation does not hold for the Levy noise 
L{t) due to the infinite covariance [2l| and L{t) can be regarded only as an external noise. The dynamical equation 
with the retarded friction, the generalised Langevin equation (GLE), is well known for the Gaussian noise. Then, 
for a given noise covariance, it allows us to determine all the fluctuations which, on the other hand, is not possible 
for non-Gaussian noises like ^{t). This process is interesting due to its jumping structure, convergent variance and 
non-trivial distributions resulting from GLE driven by ^{t). In this paper, we discuss those distributions for two 
different memory kernels: exponential and 1/t. 

The paper is organised as follows. In Sec. II, the autocorrelation function for the Ornstein-Uhlenbeck process 
with the multiplicative Levy noise is derived both for the stable and truncated distribution. Sec. Ill is devoted to 
GLE driven by that process for the Cauchy distribution: the probability density distributions are simulated and the 
resulting fluctuations are compared with general analytical predictions. A similar analysis is performed for the case 
of the power-law covariance. Results are summarised in Sec. IV. 



II. AUTOCORRELATION FUNCTION FOR THE MULTIPLICATIVE ORNSTEIN-UHLENBECK 

PROCESS 

Dynamics of a massless particle subjected to a stochastic force and a linear deterministic force is determined by 
the following Langevin equation 



(2) 



where the ^-dependent noise intensity G(^) accounts for a nonhomogeneous form of the stochastic activation. We 
assume that increments of the stochastic force, ?7(i), possess the stable and symmetric Levy distribution defined by a 
characteristic function exp(— where a (0 < a < 2) is a Levy index and G(^) = |^|~^. The system ([2]) resolves 
itself to the ordinary Ornstein-Uhlenbeck process if = and a = 2. In the multiplicative case, we must settle the 
stochastic integral interpretation that is decisive for the existence of the second moment: in the Ito interpretation the 
asymptotic distribution of ^ is the same as for rj whereas in the Stratonovich one the multiplicative factor essentially 
modifies the tail ■ The latter interpretation applies when the white noise is regarded as a limit of the correlated 
noise and the inertia is small [T2j . Then the Langevin equation ^ can be reduced to an equation with the additive 
noise by a simple change of the variable. 
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ler'sgn(e). 



K{l + 0) 

The Fokker-Planck equation in the new variable takes the form 



and its solution, after transformation to the original variable, reads 



1 + ^1.1 



1^1 



i+e 



(I, I/a), (1, 1/2) 
(1,1), (1,1/2) 



where 



a{t) 



I - exp[-7a(l + d)t] 



(3) 



(4) 



(5) 



(6) 



and the initial condition p(C,0) — (5(^) has been assumed. Asymptotics of Eq.dS]) is a power-law: p{£,,t) ~ 
1^1^^^"^"^ (1^1 1). In the limit t — > cx), the system reaches a stationary state which is characterised by the 
variance 



(^2^ ^ _^^2/(l+0)^-l-2/a(l+e)(^ _^ ^)(2a-2)/a(l+e)^-2/a(l+e)p (^-^^ T (^1 + siu [jTe) 



(7) 
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if 6* > 2/a- 1 (a < 2). 

For a — 2, a, formal expression for the correlation function can be derived by means of an expansion of the general 
Fokker-Planck equation solution into its eigenfunctions p2j |. The asymptotic behaviour appears exponential and 
the rate is given by the lowest eigenvalue. However, this conclusion may be wrong if a continuous spectrum is not 
negligible. This happens, for example, for the linear problem {9 — —1) and then the exponential is modified by an 
algebraic term |23]. The system ([2|) for a < 2 and with 6 = has a finite relaxation time but its quantification in 
terms of the covariance function is possible only after either a modification of this quantity ^23] or by introducing a 
cut-off in the distribution. In the case of the multiplicative noise, for 9 > 2/a— 1, the autocorrelation function exists 
and can be expressed by the integral 



m - mm) 



J ya6p(6,i;Ci,0)dad6- J yei6p(6,t|ei,0Mei)deid6 



where p(^) = limt^acPH, t). In terms of the transformed variables, C{t) assumes the following form 
C{t) = [K{l + 9)]^ / / \yi\^\y2\^sgii{yi)sgn{y2)p{y2,t\yi,0)p{yi)dyidy2. 



The conditional probability in Eq.® is given by 



Piy2,t\yi,0) 



H. 
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m - yie 



(l-l/a,l/a),(l/2,l/2) 
(0,1), (1/2, 1/2) 



(8) 



(9) 
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where 79 = 7(1 + ^), whereas 



\1/q 



p(yi) = — H. 
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\yi\ 



(l-l/«,l/a),(l/2,l/2) 
(0,1), (1/2, 1/2) 



(11) 



The integral ^ can be estimated in the long-time limit; details of the derivation are presented in Appendix A. The 
final expression reads 



C{t) w AKwi{l + 9)Wa-' 



■1+0 1-^^ 



1 rf^) 



V l+e) 



-/e" 



(12) 



where / is given by Eq. (jA4p . However, the expansion (lAip contains infinite terms if a < 2; in particular, the asymptotic 
expansion of the Fox function in Eq. (jA4l) produces integrand tail of the form y^^^^^^ " and then / diverges for any 
9 it a < 1. Therefore, the approximation (fT2|) is not valid for the general stable distributions. However, one can 
argue that in many systems very long jumps do not emerge and it is reasonable to introduce a truncation of the 
distribution. Such a truncation can be realised as a simple cut-off or by inserting a fast-falling tail; typical forms 
are the exponential [2^ and a power-law [26l . \2T\ . Systems involving the multiplicative Levy noise were considered 
from that point of view in Ref.[28|. Convergence to the normal distribution, expected in this case, is so slow that it 
is not observed in the numerical simulations. The cut-off at some value of \r]\ implies a finite upper integration limit 
in Ea. (|A4p and / becomes convergent. Then the autocorrelation function falls exponentially with time; the rate does 
not depend on a and rises with 9. The above result is valid also for the additive noise, 9 = 0. 

The following analysis is restricted to the Cauchy distribution of the noise rj{t) {a = 1). The integral ([9]) for the 
case without the truncation has been evaluated numerically. Inserting the conditional probability 



p(y2,t|yi,o) 



1 - cxp(-7et) 



■^le iy2 - yi exp(-70i))2 + K^{1 - expi-jot))^ / j'^ 



(13) 



to Eq.® yields the expression for C{t); it is presented in Fig.l for some values of 9. The figure reveals a stretched 
exponential shape, exp(— At''), and the parameter /3 rises monotonically from 1.040 for 9 = 1.5 to 1.081 for 9 = 5. 
Since /? is close to 1, deviation from the simple exponential emerge only for very small values of C{t) (large t) and/or 
large 9. Therefore, C{t) can be reasonable approximated by the dependence 



"-^'^ cos(7r/(l + 0))'' 



At 



(14) 



where C(0) follows from Eq.© and A is a parameter. Results for the truncated distributions, also presented in the 
figure, exhibit the fast-falling exponential tail, in agreement with Eq. (|12p . and they coincide with the stable case at 
small t. 
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FIG. 1: (Colour online) The autocorrelation function of ^{t) for a = 1 and 7 = 1, calculated from Eq.((9]), for the following 
values of 6: 1.5, 2, 3, 4, and 5 (points, from right to left). Results for the distribution truncated at lO" and 10^ are marked by 
squares and triangles, respectively. The red solid line represents the stretched exponential function and the dashed green line 
has the shape e"^"*. 



By rescaling the noise ri{t) we get rid of K and then A is completely determined by 9 and 7. Fig. 2 demonstrates that 
dependence on both parameters is simple; the expression for A(^,7) can be uniquely determined from those results: 

A = 0.80(61 + 0.31)7 = ce7. (15) 

The formula ([HI) with A from Ea. ([T5|) and (^^) from Eq.iO will be applied in Sec. Ill as the approximation to Eq.([3]) 
if time is not very large. 



III. MEMORY EFFECTS IN THE DYNAMICS 



When the dynamics proceeds in a medium of a nonhomogeneous structure, one can expect nonlinear effects and 
non-Gaussian distributions. For example, a Langevin description for the case of a Brownian particle interacting 
with a general non-Gaussian thermal bath resolves itself to a nonlinear, multiplicative Langevin equation with a 
non-Gaussian white noise and nonlinear friction term which follows from the detailed balance symmetry [29j . If the 
equilibrium state of a stochastic system results from an interplay between an internal noise and damping, the noise 
intensity and the dissipation have to be mutually related (the Einstein relation). For a correlated noise and a linear 
coupling in the thermal bath, that relation requires a retarded friction in the Langevin equation which then becomes a 
linear integro-differential equation [30I [sj- Memory effects are important also for processes involving the Levy stable 
noise. It has been demonstrated in Ref. that the external noise relaxation time modifies a slope of the power-law 
density distribution. The memory makes the friction term nonlocal in time. The fractional Langevin equation was 
introduced by Lutz 32] for the Gaussian noise which is distinguished due to the central limit theorem. For more 
general cases, e.g. in complex systems, an ordinary central limit theorem is no longer valid and the effective random 
force may assume a form different from the Gaussian even if a coupling within the thermostat is linear. GLE may 
be applied to such non-Gaussian processes [s^l but in this case higher moments cannot be expressed by the first and 
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FIG. 2: (Colour online) Dependence A(7) and X{9) (points). The straight lines mark the functions 1.857 8.nd 0.17 + 0.86 for 
the left and right panel, respectively. 

second moments. We assume that dynamics is governed by a Langevin equation with the retarded friction and driven 
by the effective random force £,{t), defined by Eq.([2|). It satisfies the second fluctuation-dissipation theorem (FDT) 
[34| to ensure a proper thermal equilibrium. Such a description is possible since FDT requires the existence of only 
first and second moments. Then we consider GLE in the form 



dv{t) 



m — = -m / K(t - T)v(T)dT + ^(t), 

at Jo 



(16) 



where v(t) is a velocity. FDT implies that the memory kernel has the same form as the noise covariance, K(t) = 
C{t)/'mT, where T is the temperature and the Boltzmann constant is set at one. The equipartition energy rule is 
satisfied: (7;(oo)^) = T/m. Applying the Laplace transformation yields the solution, 

v{t) ^ R{0)vo + m-^ [ R{t - T)^{T)dT, (17) 
Jo 

with the initial condition v(0) = vq, where the Laplace transform of the resolvent R{t) is given by the equation 



R{s) = l/[s + K{s)] 



(18) 



All the fluctuations, if they exist, are determined by the resolvent R{t). The energy equipartition rule follows from 
Ea. ([T7| in the limit t — > oo. The resolvent R{t) has an interpretation of the velocity autocorrelation function, 
Cv = {"^ovit)) = {T/m)R{t), and it determines a speed of the relaxation to the equilibrium [ssj : 



{[v(t)-R{t)v^) = -[l~R'{t)]. 

m 



(19) 



We assume that the driving noise £^{t) is given by Eq.([5]) and approximate its autocorrelation function by the 
exponential dependence (fT4| . Then a straightforward calculation yields 



R{t)^ 



'l/{2^f^) {Be 



-At 



sm 



At + V^cos- 



(A<0) 
/Ail (A>0) 



(20) 
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FIG. 3: (Colour online) The velocity distribution for 6 = 2 and 7 = 1 at t = 0.2, 0.5, 1, 2, 5, and 10 (from left to right), calculated 
from GLE, Eq.fTS)). The straight red lines mark the dependence v~'^. Inset: velocity increments for a single trajectory. 

where A = A/2 - V-A, B = A/2 + V-A and A = (^2) /(mT) - X^/A. One could expect that the probability density 
distribution, Pv{v, t), converges to the normal distribution due to the finite variance. According to Eq. ([T7|) . the velocity 
is a linear combination of the weighted values of the noise, R{t—ti)^{ti)At, where At is a constant integration step. The 
subsequent components are not independent but, since the autocorrelation C(t) falls with time, terms corresponding to 
times larger than some relaxation time of £,(t), ~ I/7, can be regarded as independent and assumptions of the central 
limit theorem are satisfied. More precisely, for a sufficiently large t the sums 5„ — J27=o ~ ~ ntr)£,{iAt + ntr), 
where n = [tr/At], are independent stochastic variables of finite variance. The variable J2n=i^^ 5'„ may converge 
with t to the normal distribution if both t/t^ is large for a non-zero R{t) and fluctuations of S'„ are small. The 
latter condition emphasises importance of higher moments of f(i). According to the Berry-Esseen theorem [s^, a 
distribution of a sum of m mutually independent variables differs from the Gaussian by 33{^'^) / {Aa^^/m), where a is 
the standard deviation, providing the third moment is finite. Therefore, convergence to the normal distribution is not 
ensured ii 9 < 2 and even for a larger 9 it may be very slow. Deviations from the Gaussian are especially pronounced 
for large values of |w| which correspond to jumps and such events are usually a result of single stochastic activations. 
In this case a distribution is similar to p(^) and tails have the form On the other hand, events that produce 

small \v\ correspond to the trajectories consisting of many small segments, similar to the ordinary Brownian motion, 
and fluctuations are small; then we may expect the Gaussian shape. 

The Monte-Carlo simulations confirm presence of the power-law tails. Time evolution of py{v,t), presented in 
Fig. 3 for the case of the infinite third moment, indicates no trace of a convergence with time to the Gaussian; the 
distribution apparently reaches a stationary state near t — 2 and the power-law dependence, dominates the 

distribution. Trajectories reveal a jumping structure typical for the Levy flights. This structure is clearly visible when 
we plot the velocity increments for the discretized integral (I17|) (Fig. 3). Note the finite jump relaxation time. 

Speed of the equilibration is governed by the parameter A in Eq. (j20p . Dependence on the parameters follows from 
Eg . (ff)) ■ ([T5|) : an estimation for large 7 yields A ^ 7~^/(^+^-'~^/mT. We conclude that the equilibration time (~ l/A) 
rises with 7 because the noise intensity declines. Fig. 4 presents equilibration of the variance for different sets of 
the parameters. The equilibrium value, (w^) = T/m, is reached at short time for small both 7 and 9 since then 
the noise intensity is large. Differences between results of the Monte Carlo calculations and Eq. p^ are due to the 
approximation of the exact C{t) by the exponential; the equilibrium value is slightly overestimated compared to the 
energy equipartition rule. 
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FIG. 4: (Colour online) Time evolution of the variance calculated from the Monte Carlo simulations (points) and from Ea. (|19p 
(lines) for T = 1 and m — 1 with the initial condition no = 0. The following cases are presented (from left to right): 1. 9 — 3 
and 7 = 0.1 (black), 2. 61 = 2 and 7 = 1 (red), 3. e = 3 and 7 = 1 (green), 4. 6* = 2 and 7 = 2 (blue), and 5. 6* = 2 and 7 = 5 
(cyan). 



The position x is given by an expression similar to Eq. pT]) but with the resolvent Rx{t) — J*R{t')dt'. The 
distribution of x, Pxix, i), is presented in Fig. 5 as a function of time; we assumed the initial condition x(0) = 0. Two 
limits, discussed above, are clearly visible if time is sufficiently long: px{x,t) assumes the Gaussian shape for \x\ < 20 
(t = 30) whereas the tail is of the form The position variance directly follows from the identity 



{x^){t) = 2 / {t-t')C,{t')dt\ 
Jo 

where the averaging is performed over the equilibrium state. A direct evaluation for the case A > yields 



{x^m 



mT 



Me] 

mT 



■t + e 



-At/2 



(A^-|^)cos(yA.) + ^(^-3g)sin(yA.) 



(21) 



(22) 



We omit the analogous expression for A < 0. In the limit of large time the variance rises linearly with time and the 
diffusion coefficient V = \imt^^{x^{t))/2t = T^X/i^) = T^ce cos(7r/(l + 0))-f^+'^^(^+'>\ The position variance as a 
function of time is presented in Fig. 6. The Monte Carlo results, obtained by integration of Eg. (fT7| . reveal a slightly 
stronger time-dependence than the linear growth predicted by Eg. ([21]) : with /3 — 1.02, 1.04 and 1.05 for 9 = 1.5, 2 
and 3, respectively. 



In many physical problems the observed covariance functions are not exponential. The power-law form of the 
memory function was discussed in connection with a frictional resistance (37j and in the hydrodynamics [ssl j. In 
particular, diffusion in the dense liguids reguires the memory function falling like i39i |. GLE for systems with 

power-law kernels, \t\~^ {f3 7^ 1), takes a form of the fractional Langevin eguation, where the damping term is 
expressed by the Riemann-Louville operator, oDt~^ [32l |: power-law kernels are present in the fractional Brownian 
motion theory J4^. Long-time correlations are observed in the complex systems that usually possess non-Gaussian 
distributions with power-law tails. A very slow falling covariance, corresponding to a 1// noise, was found in an 
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FIG. 5: (Colour online) Upper part: the position distribution for 6 — 2 and 7 = 1 at t = 0.2, 0.5, 1, 2, 5, and 30 (from left to 
right). The red dashed lines mark the dependence and the Gaussian. Lower part: the same but for the covariance (|23|) 
and the following times: 1, 2, 5, 20 and 50. 



analysis of absolute returns in the US market [4]|. The autocorrelation function of the displacement may even rise 
with time; this effect was experimentally demonstrated for the diffusion in a dusty plasma liquids for which the 



corresponding probability distributions exhibit fat tails 42|. The covariance in turn, was observed in connection 



with the noise-induced Stark broadening [43| and obtained, for a two-dimensional system, from the Navier-Stokes 
equations [13] . The presence of this form of the autocorrelation function may be related to a specific topology of the 
medium: it emerges when the trajectory has a structure of long straight-line intervals, like for the Lorentz gas [45|, 
and may be encountered in the nuclear reactions [46!|. In this paper, we solve GLE with the memory function in the 
form 1/t. More precisely, we assume 

C{t) = {l-e-^')/Lt, (23) 



where L =const> 0. 
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FIG. 6: (Colour online) Position variance obtained from integration of Ea. (|17|l with the initial condition vq = a::(0) = for 
7 = 1 and 9 = 1.5, 2, 3 (solid lines from bottom to top). Results marked by dashed red lines follow from Eq. (|21|l . 



According to the results of Sec. II, any process £,{t), given by Eq.([2]), is characterised by the exponential covariance 
and the rate is uniquely determined by A. Assuming that £^{t) is an elementary process ^x{t), we can construct a 
compound process by a superposition oi ^x{t) where the parameter A is regarded as a stochastic variable. Therefore, we 
may obtain an arbitrary, a priori assumed covariance by averaging, with a weight "tpiX), over an ensemble of trajectories 
corresponding to a fixed value of 9 and different values of 7. Since, for a given 9, C{t) = (^^)'(/;(A)e^'^*(iA ^ 
A~^/'^^+^V-'('^)e^^*c?A, the distribution of A can be evaluated for any C{t) by inversion of the Laplace transform: 

^(A) - A2/(i+«)/:-1[C(0]. (24) 
Eq. ()23|) corresponds to the following normalised distribution 

^(A) ^ -JL!l!l!^c^y(^+'^) cos ^A2/(i+«) (25) 
' 2/(1 + 6*) + 1 " 1 + 9 ^ ' 

for A € (0,L) and elsewhere. The distribution of all the elementary processes ^\{t) has the same asymptotic form, 
1^1 with the same slope of the tails since 9 is fixed in the statistical ensamble. 7, in turn, influences a relative 

intensity of the noises £,\{t). Solution of GLE is given by Eq. ([T7|) where transform of the memory kernel, K{s) = 
[ln(s/L + 1) - \n{s/L)]/L, follows from Eq.dUl). Eq.dST]) yields V = limj^oo d{x^{t))/dt = C^{t)dt - ^(0) = for 
any L. Inversion of the transform yields 

i?(t)=e-(c,cos(M) + c.sin(M))-Xy^ _ _ ^ ■ (26) 

Details of the derivation and values of the coefficients, as well as some remarks about the numerics, are presented in 
Appendix B. 

The shape of the stationary velocity distribution for the covariance ([23|) is similar to that for the exponential 
covariance case but the dependence l/|wp"^^ of the tails for 9 = 2 shifts to the relatively large \v\ and the equilibration 
time is larger. The damping parameter a, given by Ea. (|Bl[) . non-monotonically depends on L but it becomes very 
small for large L and the time needed to reach the stationary state is then extremely long. The parameter 9 does not 
influence the equilibration time but strongly modifies the distribution tail. For example, the tail assumes the shape 
for 9 = A, i.e. it falls stronger than the noise distribution. Anyway, a convergence to normal distribution is not 
observed. 
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FIG. 7: (Colour online) Position variance obtained from integration of Ea. (|17p with the initial condition vq = a::(0) = for the 
covariance (|23|l . calculated with the parameters L — 1 and 9 — 2 (points). Results marked by the red line follow from Ea. (|21|l . 
Inset: the resolvent Rx{t) for L = 0.5 (dashed line), 1 (solid line) and 2 (dotted line). 



Distribution of the position was calculated by means of the integrated resolvent Rx{t) and results are presented in 
Fig. 5. The x-dependence of Pxix, t), in particular the asymptotics, is similar to the case of the exponential covariance. 
The main difference consist in the expansion speed: the distribution for the covariance (|23p widens much slower. 
Time-dependence of the variance is given by Eg. d?!]) . where the integral in I^E^ has to be estimated numerically. 
The shape of the curve, shown in Fig. 7, reveals an apparent shape at the long time which indicate the sublinear 
behaviour. As expected from the equation R{s = 0) = 0, the system is subdiffusive and simulations agree with Eq. ((2T|) 
in the stationary limit. However, asymptotics of the variance is in fact not a power-law. According to a conjecture 
in Ref. [J?], the position variance should behave in the limit of long time like t/ln{t) the form of which has been 
interpreted by the authors as an analogy to critical exponents in a phase transition. One obtains a similar dependence 
when 'D{t) is estimated by establishing lower and upper limits of the integral [i^. The expression (x^) (t) — 0.91t/ ln(i) 
agrees with the exact result. Eg. (|2T|) . for t > 30. The time-dependent diffusion coefficient is given by the resolvent 
'D{t) ~ Rx{t) and it is also presented in Fig. 7. It appears very sensitive on L: oscillations, being strong for small L, 
vanish quickly if L is large. 

IV. SUMMARY AND CONCLUSIONS 

The overdamped Langevin equation with the quadratic potential and the multiplicative Levy stable noise describes 
a process that comprises a jumping structure of trajectories and convergent moments: variance and covariance. We 
have demonstrated that the autocorrelation function C{t) for the truncated distribution falls exponentially, in the limit 
of a long time, with the a-independent rate 7(1-1-6'). Correlations were studied in detail for the Cauchy distribution. It 
has been found that for the stable case C(t) obey the stretched-exponential form but can be reasonable approximated 
by the simple exponential. The rate has been uniquely determined as a function of the system parameters: it rises 
linearly with 7 and 9. Higher moments may also be convergent if one chooses a sufficiently large 9. The exponential 
decay for the truncated case is faster than that for the stable distribution, a conclusion that emphasizes a role of 
very long jumps in preserving the memory in the system. One may construct a stochastic process characterised by an 
arbitrary form of the covariance by a superposition of trajectories with different 7, i.e. by assuming the parameter 7 
as a stochastic variable. Moreover, one can reproduce an arbitrary slope of the distribution tail since that is governed 
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solely by the parameter 6. The above properties of the process ^ {t) suggest its applicability to problems which require 
both fat tails and long correlation time. 

If a stochastic force that obeys the above properties is balanced by the damping force, the fluctuation-dissipation 
theorem and the equipartition energy rule are satisfied; then the process obeys GLE, a fact that is well-known for 
the Gaussian case. We applied GLE to the case for which the driving force is given by S,{t). The equation predicts 
tails of both velocity and position distribution, dominated by single jumps, of the same form as the driving noise. 
The central part of the distribution, in turn, results from many small stochastic activations and for px{x,t) converges 
to the Gaussian, whereas the intermediate region assumes the fast-falling power-law. Similar distributions may be 
observed in many complex systems since they are characterised by a substantial memory and the thermal equilibration 
is accompanied by rare but spectacular events. Transport properties of the system described by GLE follow directly 
from the noise covariance; the position variance rises linearly for the exponential covariance and sublinearly for the 
covariance ^ \/t. Numerical trajectory simulations involving the process f (i) confirm that general result. Therefore, 
jumps and power-law tails of the distribution may coexist with the thermal equilibrium. 



APPENDIX A 



We derive the expression for C{t) in a limit of large t, Ea.([T^. First, we expand the conditional probability, Ea.([TUl). 
in powers of e = e"^"* to the first order. Expansion of the first term in Eq.([TU]) yields (1 — g")-i/" = 1 -|- e"/a -|- . . . . 
The Fox function is given by the series, 



H\ 



1 



\V2-Vi^\]^H[ 



("76 ) 



1/a Q 

-\y2\] + j-H{e^Q)e + .. 

Of. 



where the coefficients are dropped. The derivative involves the Fox function of the higher order pol [soj : 

(0,l),(l-l/a,l/a),(l/2,l/2) 



dH 
'd7 



yi 



2/2 - 2/ie 



rrl,2 
^3,3 



1 



12/2 - 2/ie| 



(0,1), (1/2, 1/2), (1,1) 



(Al) 



(A2) 



Next, we insert the above expansions to Eq.® and neglect terms of a higher order than e. The first component 
vanishes because the double integral can be factorised and both integrands are odd. The integral over 1/2 resolves 
itself to a Mellin transform: 



^-«/(i+.)^i,2[(a7^^^j^^^ = 1 ) = _^i/(i+«)(„^,)-i/«(i+«)£0 , (A3) 



2/2 



K 



where x stands for the Mellin transform from H^'^. Elimination of the algebraic factor in the integral over yi yields 

\ ^ a{i+s) ' a/' v-^ ^ 2 2/ 

Vi+gj ^J^ 2 1+6" 2/ 



/ = 



H. 



1,1 



dx. 



(A4) 



APPENDIX B 

We derive the expression ([26)) where the integral /J^/^ R{z)e^^dz is to be evaluated. The contour consist of a 
straight line parallel to the imaginary axis at a > 0, a large half-circle in the left half-plane and a cut along the real 
segment (— L, 0). Roots of the equation 

Lz + \n{z/L+l)-\TL{z/L) = Q (Bl) 

are of the form zi 2 — a±bi and they have to be found numerically for a given L. After a straightforward evaluation 
of the sum over residues, we obtain the first component of Eq. ipS]) where the coefficients are Ci = 2[a'^ + 2a?L + 
a{29 - l)L + a2(262 + L'^ - l) + b^{b^ + L"^ + 1)]/A, C2 = 26(2a + L) / A aivd A = \ + + h"^ + 2a^L + 2a{9 - l)L + 
b^{2 + L^) + a^(26^ + L'^ ~2). Contribution from both branches along the cut resolves itself to the integral in Ea.(l26|. 
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Trajectory numerical simulations require a value of R{t) for each integration step. Approximation of the asymptotics 
is easy to determine. For example, foi 9 = 2 and L = 1 we get R{t) = — 0.113f~^-^^^ {t > 60), a formula that coincides 
with the numerical integration up to at least t = 2000. R{t) for small t was evaluated with a step 0.001 and stored. 
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